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ABSTRACT 

The origin of crystalline grains in comets and the outer regions of protoplan- 
etary disks remains a mystery. It has been suggested that such grains form via 
annealing of amorphous precursors in the hot, inner region of a protoplanetary 
disk, where the temperatures needed for such transformations were found, and 
were then transported outward by some dynamical means. Here we develop a 
means of tracking the paths that dust grains would have taken through a diffu- 
sive protoplanetary disk and examine the types and ranges of environments that 
particles would have seen over a 10 6 year time period in the dynamic disk. We 
then combine this model with three annealing laws to examine how the dynamic 
evolution of amorphous grains would have led to their physical restructuring and 
their delivery to various regions of the disk. It is found that "sibling particles" - 
those particles that reside at the same location at a given period of time-take 
a wide range of unique and independent paths through the disk to arrive there. 
While high temperatures can persist in the disk for very long time periods, we 
find that those grains which are delivered to the cold outer regions of the disk 
are largely annealed in the first few xlO 5 yrs of disk history. This suggests that 
the crystallinity of grains in the outer disk would be determined early and re- 
main unchanged for much of disk history, in agreement with recent astronomical 
observations. 



Subject headings: astrochemistry; comets: general; meteorites, meteors, mete- 
oroids; methods: numerical; protoplanetary disks 
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Introduction 



Evidence for the large-scale redistribution of materials during the earliest stages of planet 
formation is found both in primitive objects in the Solar System as well as astronomical ob- 
servations of protoplanetary disks. The STARDUST spacecraft returned crystalline silicates 
and refractory grains which are chemically and isotopic ally similar to those materials found 



in chondritic meteor ites, pointing to a common origin (IBrownlee et al.l 120061 ; ID. et al.ll2006 



Zolensky et al.l 120061 ) . As such, it is thought that these high-temperature products formed 
in the solar nebula close to the Sun, where the temperatures needed for those objects to 
form (>~1000 K) were naturally expected, and were subsequently carried outwards by some 
dynamic process. Given that crystalline silicates are observed in the cold, outer regions of 



disks around other stars (e.g. Ivan Boekel et a 



Sargent et al.ll2009l ; ISicilia-Aguilar et al 



2009 



200 4 lApai et al. 



Oliveira et al 



2005; Watson et al. 2009 



20 111 ) -environments similar to 



that in which our comets are thought to have formed-it is thought that the same dynamical 
processes are also responsible and are natural consequences of protoplanetary disk evolution. 

A number of mechanisms have been suggested as driv ing this redistribution of materials 
in pr o toplanetary disks, including turb ulent diffusion (e.g. lGailll2001t iBockelee-Morvan et al. 



2002 



Cuzzi et allbooa Icieslal boiOah. spiral arms in gravitationally unstable disks (iBossI 



20081 ). photophoresis dMousis et al.l I2007T) . l arge-scale flows associate d with mass and an 



gular momentum t ransport (IKeller fc Gaill 12004 
acauet et al.l 2011 ). bipolar outflows and jets (e.g 



Ciesla 



2007 



Shu et al. 



20091 . l2010al : iDeschl 12007 



19961 ). and radiation pressure 



( Vinkovi?f 2009 ) . Each of these processes has been shown to carry materials from the hot, 
inner regions of a protoplanetary disk to the cooler, outer regions if appropriate conditions 
are present in a protoplanetary disk. However, it remains unclear whether these various 
conditions were met and thus which of these was primarily responsible for the large-scale 
redistribution of materials in our solar nebula. 

Regardless of the exact cause of this redistribution, the fact that solid materials get 
transported from one location to another in a disk implies that they pass through a variety 
of chemical and physical environments prior to their incorporation into a comet, asteroid, 
or planetesimal. As solids move through these environments, they would be destroyed or 
altered, the levels to which depend on the chemical and physical reactions that take place 
and the time spent in a given environment. Thus grains are not simply products of a 
singular disk environment in which they originated, but rather record the integrated path and 
physical conditions to which they were exposed over their lifetime in a protoplanetary disk. 
Understanding the chemical evolution of protoplanetary disks and the first building blocks of 
a planets requires understanding the types of paths that grains are able to take through a disk 
during the early evolution of a planetary system. As the fingerprints of such environments 
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may be found in primitive materials in comets and meteorites, this understanding may help 
us determine which one or ones of the dynamic processes identified above operated in the 
solar nebula. 

The motivation of this study is to develop a method for calculating representative par- 
ticle paths through a diffusive protoplanetary disk over long time periods (>10 5 years) and 
to apply that method to understand the origin of materials found in chondritic meteorites 
and comets. We adopt the a-viscosity model to describe the dynamic evolution of the disk, 
that is, the mass and angular momentum transport and the diffusion of species within it. 
The a- viscosity model is used because it provides a means of describing the evolution of the 
protoplanetary disk and the materials within it over timescales of millions of years, which 
is necessary when examining the data provided by meteorites and primitive materials in the 
Solar System. This model assumes that the evolution of the disk is tied to some turbulent 
viscosity, u=ac s H, where c s is the local speed of sound, H is the local scale-height, and 
a is some parameter (a <1) which quantifies the strength of the turbulence of the disk. 
The relation of radia l diffusivity to this viscosity remains the subject of ongoing work, (e.g. 



Johansen et al.l 120061 ; iPavlyuchenkov fc Dullemondl 120071 ) , and we will begin by making the 
standard assumption that the diffusivity of gas in the disk is D g ~ is, though the techniques 
described here do not require this assumption. 

In the next section, the particle-tracking model is described and validated against other 
calculations for the radial transport of materials in a protoplanetary disk. In Section 3 we 
describe the types of paths particles take through a turbulent, viscous protoplanetary disk. 
In Section 4 we use these results to explore how the motions of dust grains in a protoplanetary 
disk affect their physical structure by feeding information about the environments they saw 
into various kinetic models for the annealing of silicates. In Section 5 we summarize the 
paper and present our conclusions. 



2. Model 



The methods used here build off of that described by ICieslal (j2010bl ) who developed an 
approach to calculat e representati ve vertical paths of particles in diffusive protoplanetary 
disks. In that study, ICieslal (j2010bl ) too k the advection-diffus i on equation that describes th e 
vertical evolution of a species in a disk (IDubrulle et al.lll995t iFromang fc Papaloizoull2006l ): 



dPi_d_l n H£ 

dt dz Pg dz 



q \PiVgra 
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where pi is the spatial density of the material of interest, p g is the gas density, D the diffusion 
coefficient, and v grav the vertical velocity of the species which arises due to gravitational 
settling and calculated the vertical motions of an individual particle using the equation: 



Zi-l + V e ff5t + R 



T D iz') St 



(2) 



where and Zi described the position of the particle before and after a timestep of duration 
St. Eq (E} came about from analyzing the Fokker-Planck equation as well as by considering 
the moments of Eq ([I]). The second term on the right of Eq. ([2]) represents the advective 
motions of the particle, with the velocity not simply being the gravitational settling velocity, 
but includes effects that arise from the spatial variations in the gas and diffusion coefficient. 
Thus, the effective velocity, f e //, is the sum of three terms: 



V eff = V grav + V gas + V F P 



(3) 



with v grav again being the vertical settling velocity set by balancing the vertical component 
of the gravitational force from the central star with the resistive drag force as the solid moves 
through the gaseous fluid 
of materials depends on the ratio Pi/p g and 



v gas is a term that arises due to the fact that the diffusive flux 



varies with z (v. 



gas~ 



. D dp g ■ 

pg dz • 



and vpp is a 



term that accounts for possible variations in the diffusi on coe f ficient D (vfp=^)- Detailed 
derivations and explanations of each term are given in Icieslal f 2010b ). 



The third term on the right of Eq. (j2J) represents the displacement that occurs as a 
result of the diffusive motions in the disk. This term is evaluated by finding a random 
number, R, from a distribution with variance £. The diffusivity is evaluated at a location 



Z =Zi 



-i+§7j?<^ to account for strong gradients in the diffusivity. 



Here we want to develop a similar means of calcu lating the radi al motions of materials 
in a diffusive protoplanetary disk as that described in ICieslal (l2010bl ). The one- dimensional 
radial advection-diffusion equation that describes the behavior of a species, i, in a pro t oplan- 
etary disk is written (e.g. iGaill l200ll ; iBockelee-Morvan et al.ll2002l ; ICuzzi et al.ll2003t ICiesla 
2010ah : 

as, i d ( „„ d m\\ is, „, ... 

(4) 



rS.D - 

dt r dr V dr V S 



^- {rv r T,i 

r or 



where Sj is the surface density of the species of interest, S is the surface density of the disk, 
D is again the diffusivity, and v r is the radial motions of the materials due to large-scale 
flows and gas drag. The form of this equation is slightly different from that of Eq ([!]) due to 
the cylindrical coordinate system, and as such, leads to different, and less useful, relations 
between the moments of the equations. As a result, this form of the radial transport equation 
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does not immediately allow us to write an equation similar to Eq. 
trajectory of a particle in a diffusive disk. 



ff2T) to describe the radial 



A solution to this problem is to instead adopt a Cartesian coordinate system and track 
the motions of the particles in the x-y plane, where the radial distance from the origin is 
given by r 2 =x 2 +y 2 . This allows us to write the advection-diffusion equation as: 



dt 



d_ 

dx 



EL> 



d_ 

dx 



d_ 

dx 



(v x Ei) + |- ( ED 



d 



5* 

E 



d_ 

dy 



VyYji 



(5) 



dy \ dy 

where v x and v y represent the x and y components of v r . The motions in each direction are 
now described using equations that have the same form as the vertical diffusion equation 
(Eq (pQ)). A similar approach is taken in other st udies of transport , such as those examining 



the dynamics of pollutants in groundwater (e.g. iKinzelbachl Il990h . and has been shown to 



accurately describe the motions of materials of interest. 



Following ICieslal feOlObl ). the motions of an individual diffusing molecule or particle 
would be described by: 



X j 



-l + v e f f , x 8t + R 



and 



Vi = Vi-i + v e ff t y8t + R 



-D ix') St 
r 



-D{y') St 



(6) 



(7) 



where each term in the above equations represents the same physical effects as described 
for Eq (j2J). Here v e ff iX and v e ff >y have multiple components as described above. That is, 
+ vpp^ x , meaning effective velocity in the x-direction would be given by 



J eff,x z 



~-V x + V 



gas,x 



the sum of three different effects: the velocity of the materials in the x-direction due to 
large-scale flows and gas drag (which normally result in purely radial motions): 



v x [x i -. 1 ,y i - 1 ) = v r 



Xi-i 



i-l 



where rj_i is the radial position of the particle at the beginning of the timestep, v r is the radial 
advective velocity evaluated at rj_i, and is its position along the contribution 
due to the gradient in the surface density: 

D dE 



Vgas,x lj Vi— l) 



(9) 



E dr Ti_i 

where D, E, and the derivative of the surface density are evaluated at the point yi-%)] 
and a contribution due to the gradient in the diffusivity: 

dD Xi-i 



vfp,x (Xi-i,yi-i 



dr rj_i 



(10) 
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The dependency in the y direction would be similar as above. Note that we are ignoring 
orbital motions here, as we are assuming axial symmetry-the only concern is variations with 
distance from the origin. Were we concerned with azimuthal variations, orbital velocities 
would have to be added to the expressions given above. 

To validate this model, we compare the results of this particle-tracking model to those 
that solve the radial advection-diffusion equation Eq. (4) using a finite-volume method. The 
so called "dye-tracking" approach, which tracks how the abundance of a species changes as 



inves t igated in previous s tudies (e.g . iGail 



2003; Keller fc Gail 2004; Boss 2008 



2001 



Ciesla 



2007 



2009 



transport has been 


2002 




Cuzzi et al. 



2010al . among many others). The 



dye-tracking m odel calculations in this study were performed using the model described in 
Cieslal fcoiOah - 



The first validation test compared the predictions of the two models for the case of 
transport in a steady-state protoplanetary disk around a solar-mass star, whose physical 
structure was given by: 



E (r) = 2000 



1AU 



T(r) 



280 



1AU 



/ 2 

g/cm 



K 



(11) 
(12) 



As we have assumed a standard a-disk in steady-state, the advective velocity everywhere 
was given by v r —-3u/2r, where v is the local viscosity. A value of a=10~ 3 was assumed, and 
all particles were taken to be l//m in radius. 

We tracked the transport of the trace species in the disk using both types of transport 
models described above. The results are shown in Figure 1. The trace material in the 
dye-tracking model was released in the grid that spanned ~5.0 to ~5.4 AU (centered at 5.2 
AU), such that Ej=l g/cm 2 at the beginning of the simulation. We then calculated how Ej 
evolved with time and location. The results for this model at t=10 4 , 10 5 , and 10 6 years are 
shown as the solid lines. 

In the case of the particle-tracking model, iV=10,000 particles were released at 5.2 AU at 
the beginning of the simulation. Each particle thus represented ylE;/iV=9.36x 10 22 g of trace 
materials, where A represents the area of the annulus spanning from ~5.0 to ~5.4 AU. In 
calculating the evolution of the trace particles, Equations (6) and (7) were used, with the local 
values of the gas determining the appropriate velocities and diffusive step. The timestep for 
each particle was assumed to be 5t=Q.l/Q,K, where Qk is the local orbital angular velocity at 
the location of the particle at the begi nning of the time step. This i s roughly the correlation 
time expected for nebular turbulence (jFromang fc Papaloizoull2006l ). Particles were released 



-7- 



at t=0 and tracked for a period of 10 6 years, with their locations recorded at 10 4 and 10 5 
years. When the timestep would have pushed the total simulation time for a given particle 
beyond the time of interest, the timestep was was reduced to allow the locations to be 
recorded. In calculating the surface density of the trace species, the number of particles 
located between different grid points used in the dye-tracking calculations were determined, 
and the total mass represented by these particles then divided by the effective area of the 
annulus in which they were located. These surface densities distributions are shown as 
dashed lines in Figure 1. 

We also considered transport in a viscously evolving disk, thus relaxing the assumption 
of a steady-state disk. For this case we started with a O.1M disk with an initial surface 
density of 

E(r) = 14250 (iXu) _1 S/ Cm2 ( 13 ) 

out to 10 AU. Again, the standard a-disk model was used, with the model of ICiesla 
used to calculate the temporal evolution of the disk along with the transport of materials 
within it. Trace materials were tracked in the same manner and with the same initial 
conditions as described in the steady-state case. These results are shown as the solid lines 
in Figure 2. 

The particle-tracking calculations were performed as described above, however, the 
changing dynamic environment within the disk had to be accounted for. The properties 
of the disk (surface density, temperature, etc.) were recorded at every grid point in the disk 
model at 100 year intervals in the case of the evolving disk. At each timestep in the particle 
tracking model, the properties of the gas and dynamics of the disk at the location of interest 
were found by linear interpolation from this spatial-temporal map of disk evolution. The 
results of the particle-tracking model in the evolving disk are again shown as the dashed 
lines in Figure 2. 

As can be seen, there is strong agreement between the two different transport models 
in cases considered here. One immediately notes that the inferred surface densities for the 
particle-tracking model at 10 6 years are not as smooth as the earlier times within each 
case. This is due to trace materials (particles) being lost from the disk over time as they 
are accreted by the central star. This results in a smaller number of particles surviving at 
these later stages. These smaller number of surviving particles are thus not sufficient to 
provide a detailed description of the collective behavior of all particles in the disk, though 
the remaining particles still reproduce the general distribution quite well. This agreement 
demonstrates that the particle-tracking model provides a good representation of the motions 
of particles that are subjected to the diffusive and advective motions expected in the model 
disks described here. 
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3. Particle Paths and Environments Seen 

The advantage of the particle-tracking model over traditional dye-tracking models is 
that it allows one to investigate the particular paths that materials would take over their 
lifetimes in a protoplanetary disk. This is demonstrated in Figure 3, which shows the paths 
of 3 particles from the simulation of the evolving disk described above. These particles 
originated at the same location in the disk (5.2 AU), and each finished the simulation (after 
10 6 years) located within <0.01 AU of one another at ~2.56 AU. As can be seen, however, 
each particle took very different paths to reach that final location, with each moving as far 
as 13, 22, and 27 AU from the central star over their lifetimes. Given that the particles 
migrate through different regions of the disk during their 1 million year disk lifetime, they 
will see different thermal, pressure, and irradiation environments from one another. This 
would allow each grain to be processed differently within the disk prior to meeting up again, 
where they could potentially accrete into a common parent body. 

Figure 4 shows the different paths of the particles during the first 10 5 years of the 
simulation, both in terms of their distances from the star, as well as the ambient temperatures 
that they would see. While two particles are pushed outwards fairly rapidly, moving from 
~5 to 15 and 20 AU in this 10 5 year period, the third particle, whose path is indicated 
by the dark, black line in Figure 4, moves around the inner region (1-8 AU) of the disk 
over this same time period. As such, the two outwardly drifting particles move to cooler 
regions of the disk and generally never see temperatures above ~150 K after 40,000 years of 
model time. The third particle, however, spends the first 10 5 years at temperatures above 
200 K, with a number of excursions into regions of the disk with temperatures above 600 
K, including seeing temperatures above 1000 K multiple times in the first 10 4 years. Such 
temperatures are high enough that grains could have been subjected to physical processes 
such as annealing, or undergone chemical alteration through rapid reactions with the gas- 
processes that the other two grains would have escaped because they never saw temperatures 
this high. 

While these concepts are not necessarily new, the paths of the three particles shown in 
Figures 3 and 4 illustrate the importance of understanding what paths particles took over 
their lifetimes in a protoplanetary disk-it is not enough to assign the chemical compositions 
or physical properties of materials based on their initial or final locations in the disk. Figure 
5 shows the range of locations seen by the 3803 surviving (out of 10 4 modeled) particles that 
reside outside of 1 AU after 10 6 years of evolution in the protoplanetary disk (and did not 
migrate within 0.1 AU of the star at any point). Plotted are the maximum and minimum 
radial distances from the star for each particle as a function of their final location in the 
disk for both the entire disk (left panel) and for those that finish the simulation within 5 



-9 - 



AU of the star (right panel). For any given location, the materials present at the end of 
the simulation have seen a range of distances from the star, and thus, physical environments 
within the disk. For example, in the asteroid belt region, roughly the 2-4 AU region of the 
disk, we see some grains have moved inwards as far as 0.2 AU in the disk, while some have 
never gotten closer than their final location of 2-4 AU. These same grains typically migrated 
outwards beyond 10 AU, with some seeing distances as far out as 40 AU before migrating 
inwards again. 

The various distances seen by the particles at a given location would suggest different 
thermal histories for those grains. Figure 6 shows the highest temperature each particle saw 
in the disk plotted as a function of its final radial location. Those particles that remain 
in the inner disk (< 5 AU) have seen a wide range of peak temperatures, from ~620 K, 
corresponding to the temperature of the disk at their point of origin at the start of the 
simulation, to ~ 2000 K, which is achieved much closer to the star (most grains would have 
been vaporized at such temperatures, but we ignore this effect for now-we return to this 
issue below). This same range of temperatures is seen in those particles that finish the 
simulation inside of 50 AU, though there is a trend that the further out in the disk that a 
particle is located at the end of the simulation, the lower the range of peak temperatures seen. 
Specifically, at greater distances from the central star, the lower the maximum temperature 
seen by the particles present there. This is due to the fact that the particles that migrate 
furthest outwards in the disk do so by having short residence times in the inner disk, and 
thus are able to ride the viscous expansion of the disk (residing in the region of the disk 
where the net motions are outward) to very large distances from the star. As a result, grains 
that start at a given location and then wind up in the inner disk will have been exposed 
to a wider range of peak temperatures, whereas those grains in the outer disk rarely saw 
temperatures in excess of the temperatures present in the disk where they originated. 

It is important to remember that the results shown here consider only the dynamics of 
those dust particles released at a given location in the disk. This point is selected to reside 
far enough out in the disk that temperatures are too low for any significant annealing to 
take place and thus the grains are assumed to begin amorphous, which is important for the 
discussion below. While the specifics of the dynamical evolution of particles released at other 
locations are expected to differ from that here, the general trends should remain the same. 

Figures 7 and 8 show the radial ranges and peak temperatures seen of particles as Figures 
5 and 6, but for a disk in which the turbulent parameter, a, was set to 10~ 4 (with the same 
initial mass distribution). Particles were released at the same location in the disk as the case 
just discussed, 5.2 AU, which is at a lower temperature than above, ~250 K due to the lower 
levels of viscous dissipation. Because the level of turbulence is lower, the random motions 
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of the particles in a given time interval are smaller than the preceding case. Thus the range 
of behavior for particles that wind up at a given location in the disk is smaller than in the 
previous case. That is, those particles that finish the simulation between 2 and 4 AU rarely 
migrated outwards of 10 AU from the central star, while still exhibiting a similar range of 
minimum radii as above. 

These results illustrate that each particle follows its own, independent path within a 
diffusive protoplanetary disks, and thus sees a unique integrated set of protoplanetary disk 
environments. As such, those particles that are found in close proximity to one another at a 
given time would have seen very different conditions within the disk, which would translate 
to different chemical and physical evolutionary histories. The particular ranges of conditions 
seen by a collection of particles will depend on the level of turbulence (diffusion) in the disk. 
The differences in disk conditions seen by "sibling" particles would not be revealed in the 
dye-tracking models traditionally used in transport studies, and can only be quantified by 
following the representative paths that particles take through the disk as done here. 



4. Application to Annealing of Amorphous Dust 



Outward transport from the hot, inner regions of the solar nebula has been a focus of re- 
cent work as it may offer a way of explainin g the presence of c rystalline silicates in comets and 



the outer regions of protoplanetary disks (|Nuth et a 



2008; Bockelee-Morvan et al. 2002; Gail 2001 



Ciesla 



2000; 



2007 



Nuth fc Johnson 2006; Wooden 



20091 ). That is, silicates that were 



present in the parent cloud core from which our solar system formed are expected to have 
been predominately amor phous, based on the high abundance of such grains in the diffuse 
ISM (IKemper et al.l 12004 ). Such grains would remain amorphous once incorporated into 
the sola r nebula, unless exposed to high enough temperatures that they could undergo an- 



nealing (IFabian et a 



Roskosz et al. 



2009 



2000; lHallenbeck et al.ll2000l ; iDiouadi et al.l 120051 : Murata et al.l 120071 : 



20111 ). The level to which such grains were transformed depends on the 
temperatures that the solid particles were exposed to within the solar nebula and how long 
they spent at the different temperatures. 

Previous studies for the formation and redistribution of crystalline grains in the so- 
lar nebula have generally accounted for annealing by assuming precursors become crystalline 
upon entering a region of the disk that was above some critical temperature, T a . This temper- 



small that it could be considered ins tantaneous (e.g 



Gail 


2001; 


Bockelee-Morvan et al. 


2002 



Dullemond et al.l 120061 ; ICieslal 120071 ). However, annealing could be a drawn out process, oc- 
curring over an extended period at lower temperatures. Further, it is unclear what the extent 
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of the annealing would be for each grain in these studies-that is, how often are there grains 
that undergo partial annealing, but do not become fully crystalline? These issues can only 
be addressed by tracking the exact temperature-time history of grains in a protoplanetary 
disk and applying the models for annealing rates that have been derived from laboratory 
studies. 

While annealing is recognized to require high temperatures, the exact temperatures 
and rates at which this transformation occurs depend on the structure and composition of 
the amorphous pre c ursor. Experimenta l studies on Mg-rich olivine and chondritic minerals 
(IFabian et al.l I200CH : iDjouadi et al.l 120051 ; iMurata et al.l 120071 ) have shown that the an nealing; 
process is well described by the Johnson-Mehl-Avrami (JMA) equation (jBurkdll965l ): 



X 



exp 



(14) 



where X is the fraction of displaced atoms that are transformed (annealed) in a given time, 
t, and n is an exponent that depends on the crystallization kinetics. The characteristic time 
constant of transformation, r, is dependent on the temperature and follows an Arrhenius 
relation: 

1 — ve"^ (15) 



where v is a characteristic vibrational frequency and E a is the activation energy. While ex- 
perimental results vary somewhat, typical values of the p arameters for amor phous Mg 2 Si04 
smokes are are n=l, z/=2x 10 13 s"\ and £ a /A;=39,100 K (IFabian et al.lhoooh. with values of 
z /=2.5 xl0 26 s -1 and £' f ,,//c=7 0,000 K estimated fo r MgSiOs glass particles (IBockelee-Morvan et al 



( 2002 ) based on analysis of Fabian et al. ( 200ol ) experiments). Importantly, Djouadi et al. 



(120051 ) demonstrated that the annealing of the amorphous silicates is independent of the 
amorphization history of the grain-that is, the annealing of the grain will proceed in the 
same manner regardless of whether the grain is completely amorphous from its lifetime in the 
parent molecular cloud or if it was cycled in and out of warm environments, partly annealing 
over multiple intervals throughout. 



Hallenbeck et al.l (120001 ) offered an alternative way of tracking the change in the struc- 



ture of Mg-silicates by observing the evolution of their spectra over t ime as they were ex - 



posed to different temperatures over different periods of time. L ike IFabian et al 



I 

they used Mg-smokes as their starting point, though the smokes in lHallenbeck et al.l (120001 ) 
were formed from v apor phase condensates that were quenched at low p ressures, whereas 



(120001 

(hoot 



Fabian et al.l (120001 ) smokes were formed at slightly higher pressures. The lHallenbeck et al. 



(boooh smokes have been referred to as "chaotic" as they require some warming and annea l- 
ing; to occur just t o get to structures similar to those smokes studied in IFabian et al.l (120001 ) . 
Hallenbeck et al.l ( 120001 ) noticed that the spectra of the annealing grains evolved in distinct 
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stages, rather than continuously. Most importantly, they described the spectral stall, where 
little change occurred in the 10 /im region of the spectra for grains heated to just above 
1000 K, but rapid evol u tion w as seen when slightly higher temperatures were reached. As 
such, lHallenbeck et al.l (120001 ) defined a two-stage annealing process that they quantified 
by defining the Silicate Evolution Index (SEI). Amorphous grains would begin with SEI=0 
and evolve with a given rate up until the stall point, or SEI=1. Beyond that, the rate of 
evolution would change, with the spectra continuing to evolve towards that of crystalline 
silicates, which is quantified with the SEI increasing up to a value of 102. 

A comparison of the three annealing rates considered here is shown in Figure 9, where 
the time constant of transformation is plotted as a function of temperature. While other 
annealing rates have been derived, they tend to be similar to or fall within the range of those 
shown here. To illustrate the importance of u nderstanding t he dis k environments seen by a 
given particle, consider the annealing rate of iFabian et al.l ( 120001 ) (solid line). Amorphous 
grains would take ~1 year to anneal at a temperature of ~800 K, and such a temperature 
may be used as an estimate for T a , the critical temperature for annealing, as this timescale 
is comparable to a typ ical dynamical timescale in a disk (e.g. this was the criterion in 
Dullemond et al.l 120061 ) . However, amorphous grains brought to a temperature of 650 K 



would also anneal if maintained at this temperature for a period of ~10 5 years. In order to 
fully evaluate how annealing proceeds in a protoplanetary disk, it is necessary to quantify 
how long particles may reside in these lower temperature environments where annealing may 
still occur on ti mescales that are less than the typical lifetime of a protoplanetary disk (of 
order 10 6 years; lHaisch Jr. et al.ll200ll ). Note we are ignoring any mineralogical changes that 
may occur, though thi s may be important i n dete rmining; the range of temperatures at which 
grains were annealed (IRoskosz et al.ll2009l . l201ll ). 



We have used the calculations in the validation cases described above to evaluate how 
amorphous precursors, starting at an initial temperature of ~620 K (for the a=10~ 3 case), 
would be physically altered as a result of their motions through disk. Each grain was assumed 
to be perfectly amorphous (Xq=0 or SEIo=0) at the beginning of the simulation. The extent 
of annealing for a grain was determined by integrating over its individual path in the disk, 
tracking the temperatures (taken as the midplane temperature) each grain saw and the time 
spent at e ach temperature. Ea ch of the annealing /cryst allization rates identified abov e, the 
MgSi0 4 of IFabian et~aJ (boooh smokes, MgSiQ 3 g lass of Ifjockelee-Morvan et all J2002h . and 
the Mg-silicate smokes of lHallenbeck et al.l ( 120001 ) were considered. The results are shown in 
Figures 10-12, which plot the crystallinity fraction achieved by each grain as a function of the 
highest temperature the particle was exposed to over its nebular lifetime (left panels of each 
figure). Also plotted in each figure are the crystallinity fraction (fraction of the displaced 
atoms that have been transformed) for each grain versus the final location of the particle at 
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the end of the simulation (set as 10 6 years). 

For each of the annealing rates, the grain distributions form a near step function in that 
grains are completely amorphous (fraction annealed or SEI of zero) if peak temperatures 
were below a certain value, transitioning to completely crystalline (crystalline fraction of 
1 or SEI of 102) over a small temperature range, of order tens of Kelvin. This is due to 
the exponential dependence of the annealing laws, which cause the annealing timescales to 
rapidly decrease with increasing temperature. The temperatures at which t he transitions 
from amorphous to crystalline t ake p lace a re ~ 700, 850, and 1000 K for the Fabian et al. 



(120001 ). IBockelee-Morvan et al.l ( 120021 ). and lHallenbeck et al.l (120001 ) annealing laws, respec- 
tively. Thus in models where kinetics are not explicitly accounted for because the paths 
of the particles are not tracked, such as the dye-tracking models of previous studies, these 
temperatures represent good approximations for the temperature at which the amorphous- 
to-crystalline transition is made. 

The step function behavior of annealing implies that there will be an almost binary 
nature of grains in the protoplanetary disk: they are generally either 0% or 100% crystalline, 
with a small fraction ac hieving intermediate values. The fraction achieving intermediate 
values is greatest in the iFabian et al.l (120001 ) law, as displaced atoms are more mobilized 
at lower temperatures. In examining grains from chondritic meteorites, interplanetary dust 
particles (IDPs), and cometary grains, there are very few that seem to suggest that they 
were only partially a nnealed (S. Messenger, person al c ommunication). This su ggests that 
the annealing laws of IBockelee-Morvan et al.l (120021 ) or lHallenbeck et al.l ( 120001 ) may more 
accurately describe the behavior of silicates in our solar system. 

This general behavior is true regardless of the size of the particles considered. Figure 13 
shows the level of annealing achieved for grains 10 //m, 100 //m, and 1 mm in radius plotted 
against their final location in the disk after 10 6 years of evolution. The dynamics of the 
various sized particles differ primarily due to larger particles being dragged inwards due to 



the intera ctions with the gas ( Weidenschillinglll977l ). with the lower diffusivites of the larger 



particles (jYoudin &: Lithwickl 120071 ) playing a secondary role in these cases. These effects 
can be seen as the 1 mm particles are all found inside of ~30 AU in the disk at the end of 
the simulation, while the outermost distance for each other particle increases with smaller 
size. These same effects also lead to a decrease in the survival frequency of the particles 
with increasing size. Despite these differences, however, the fraction of grains that survive 
that are completely annealed is constant within the disk. Further, there still remains only 
a small number of particles that undergo intermediate levels of annealing, again implying 
grains will largely be either perfectly amorphous or perfectly crystalline. 



It is worth noting again that vaporization of dust grains is not explicitly modeled here. 
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Olivine and pyroxene grains, the dominant silicates found in primitive bodies in our Solar 
System, are expected to vaporize at temperatures of ~1300-1500 K within the solar neb- 



ula, with the exact v a lue depending o n the pressure of the gas (jEbel fc Grossman! |2000| ; 
Davis fc Richterl 120051 ; iGrossmanl |2010| ). In reality, the grains that reached these temper- 
atures would vaporize, and could only be preserved in the disk if the resulting vapor was 
carried outward again to recondense. Indeed, a number of those grains in the model devel- 
oped here saw temperatures of T >1500 K, yet remain in the disk after 10 6 years, implying 
this would have occurred. Upon condensation, the silicates may have been amorphous rather 
than crystalline, though they are treated as being crystalline in the model. This is justified 
as the condensation temperature for the silicates would largely exceed 1000 K, so that even 
if amorphous grains condense, they would quickly anneal as shown here. If it were only va- 
porization and recondensation that allowed crystalline silicates to form, the distributions of 
crystalline grains shown here would be similar, though with a smaller fraction of crystalline 
materials. Plots of the crystallinity versus peak temperature would continue to remain as a 
step function, with the transition of amorphous to crystalline occurring right at the temper- 
ature where vaporization occurred (generally higher than the transition temperatures in the 
annealing laws here). 

In addition to modeling the kinetics of the annealing process, we can also examine when 
the crystalline grains in the disk were annealed. Figure 14 plots the time of crystallization 
of the grains, defined as the time in the model when a grain reached 100% annelaing or 
SEI=102, as a function of distance from the central star. As can be seen, silicates are 
annealed throughout the entire 1 million years of model time considered here -as long as 
there is a region in the disk that is sufficiently hot and particles are able to reach it, then 
annealing can occur, regardless of the annealing law used. A trend can be seen in that those 
grains in the outer part of the disk are annealed very early on in disk evolution, whereas those 
grains that are annealed at later times are found closer to the star. This is again understood 
in terms of the dynamical evolution of the disk-those grains that migrate furthest outwards 
in the disk do so early on by taking advantage of the viscous expansion of the disk. At 
later times, outward transport becomes more difficult, as the gas has a net inward flow for 
a larger radial expanse of the disk. Thus any grains that are exposed to high temperatures 
in the inner region of the disk at later times are only able to diffuse outwards a relatively 
small distance before the inward flows dominate their motions. This same effect is what 



Cieslal (I2010af) argued was responsible for the narrow age distribution of CAIs in chondritic 
meteorites ( jThrane et al. 2006 ; Jacobsen et al. 2008 ). 



This implies that the crystalline grains observed in comets t o day s hould have all been 
annealed early in the formation of the solar system. iNuth et al.l ( 120001 ) proposed that the 
ages of comets could be inferred from the relative proportion of crystalline and amorphous 
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silicates: ancient comets would contain largely amorphous grains, while those that formed 
later could incorporate more and more annealed grains brought outwards from the inner 
solar system. Here we see that while grains may be delivered to the outer disk over an 
extended period of time, it is most likely to occur early in disk evolution. That is, those 
grains that are delivered to the outer disk likely annealed during a short period of time (<10 5 
yrs) very early in the evolution of the disk, because those grains that migrated to the high 
temperature region of the disk early on could anneal and be pushed out to where the comets 
formed by the viscous expansion of the disk. Those grains that originated in the cool, outer 
regions of the disk would only anneal once they migrated to the hot, inner disk and could 
only be incorporated into comets if they were transported to the outer disk again. This 
process becomes less and less efficient over time as discussed above-the dynamical evolution 
of the disk limits how much material can make it to the outer disk as it evolves. 

This behavior may result in a difference in the amount of crystalline silicates that are 
incorporated into Oort cloud and Kuiper Belt comets. That is, Oort cloud comets are 
thought to have formed in the vicinity of Jupiter and were subsequently scattered outwards 
through gravitational interactions with the giant planets. The Kuiper Belt comets, on the 
other hand, largely formed outside the orbits of the giant planets. These comets thus may 
contain different populations of annealed grains, with those in the Oort cloud able to accrete 
a portion of the grains that were annealed later in disk history, whereas the Kuiper Belt 
comets would only be able to accrete those grains that were annealed earlier. 

This result serves as a prediction of the annealing model for the source of crystalline 
grains in the outer disk-in the context of the model described here, all such grains should 
have been exposed to high temperatures very early in the evolution of the disk. If high tem- 
perature minerals instead exhibit a wide range of ages, including some that were exposed 
to such temperatures a few xlO 5 years into the evolution of the solar nebula, an alternative 
source of crystalline grains ma y be needed, such as th e possible annealing of grains in tran- 



sient heating events like shocks (IHarker fc Deschll2002l ) . While dating the timing of annealing 



of such grains may b e difficult, it is worth considering the recent observations reported by 



Oliveira et al.l (120111 ). who compared the dust mineralogy in a number of protoplanetary 
disks in clusters of different ages. It was found that all disks have roughly the same silicate 
crystallinity fraction, independent of age. They used this observation to suggest that anneal- 
ing of silicates would have been limited to very early in disk evolution <1 Myr, and then the 
crystallinity was "frozen in" for the rest o f disk evolution, remaining roughly constant with 



time. A similar conclusion was reached by ISicilia-Aguilar et al.l (120091 ). based on their study 



of the 7] Chamaeleontis cluster. This is consistent with the picture developed here, where 
the crystalline grains form early on in the high temperatures in the inner disk, and then 
migrate outwards to be preserved in the disk. The crystalline fraction would not increase 
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significantly as it would become more difficult for newly annealed grains to be preserved in 
the disk at later times as the dynamical environment in the disk would continuously push 
these grains to the inner disk. Indeed, the results shown here suggest that the youngest 
crystalline grains are found only in the inner disk, allow ing the outer d i sk to remain at a 
roughly uniform crystalline fraction, which is the region lOliveira et al.l ( 120111 ) was largely 
probing. 



5. Discussion and Summary 

Evidence for mixing and transport is abundant in the Solar System. Chondritic mete- 
orites contain an array of materials that formed in different environments yet are mixed on 
the scale of centimeters. These same types of materials were also found in the Comet Wild 
2 samples returned by Stardust, pointing to a large-scale exchange of materials between the 
inner solar nebula, that is the region in which terrestrial bodies formed, and the outer nebula 
where icy bodies and giant planets formed. That small solids underwent such large-scale ex- 
cursions in the early Solar System implies that they would have passed through a variety of 
nebular environments within which the solids would have had their chemical, isotopic, and 
physical properties defined or altered prior to their incorporation into a meteorite parent 
body, comet, or planetesimal. 

Here we have provided a means for calculating representative paths that solid particles 
would have taken within a diffusive protoplanetary disk. This allows us to identify the 
specific environments seen by a given dust particle and quantify how long it resided in 
each. Such information is critical to understanding the chemical and physical evolution 
of a protoplanetary disk, and as such, the method developed here provides a new way of 
quantitatively studying the origins of primitive materials. We find that within a turbulent, 
viscous protoplanetary disk, diffusion causes grains to follow a wide range of independent 
paths. 

We have illustrated how this method can be used by applying it to under stand the an- 
nealing of amorphous silicates in a protoplanetary disk. T hree annealing rates ( jFabian et al. 



20001 ; iBockelee-Morvan et al.ll2002l ; lHallenbeck et al.ll2000l ) were used to explore how the dy- 



namical evolution of individual grains impacted their physical properties. We found that 
amorphous grains largely remain completely amorphous or become perfectly crystalline, and 
that rarely do grains fall in between. This is true for each of the different annealing laws 
investigated in detail here, suggesting this is a robust result. Those amorphous grains that 
become annealed generally undergo this transformation early in the evolution of the disk. 
This is particularly true for those that wind up far from the star and could possibly be 
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accreted into comets, as they require the viscous expansion of the disk to aid their outward 
transport. That is, crystalline silicates in the outer disk likely formed over a relatively small 
window in disk history if they formed via annealing in the hot, inner solar nebula. If their 
crystallization ages could be determined, this would serve as a test for this model for their 
formation. If the crystallization ages spanned a wide range of time or were found to come 
from a late period in the history of the disk, this would suggest transient heating events 
were likely a factor in their origin. We have considered a number of other disk structures, 
evolutionary parameters, and starting locations not shown here and find that these results 
are robust. T his may explain the r eported invariance of silicate crystallinity with disk age 
as reported by lOliveira et al.l ( 1201 lh . 



The disk model used here used the standard a-viscosity prescription to describe the 
physical evolution of the disk as it underwent mass and angular momentum transport. The 
disk model thus underwent relatively "smooth" evolution, with the mass accretion rate of 
material from the disk to the central star monotonically decreasing with time. Observations 
of protoplanetary disks suggest that accretion, particularly early on, was more episodic, with 
the rate at which material falls onto the star varying by orders of magnitude on tim escales 
short compared to the lifetime of the protoplanetary disk (IHartmann fc Kenyonlll996l ). Such 
outbursts have recently been modeled in the context of an a-viscosity disk, where a var- 
i ed with location and time, as its value was depe ndent on the loc al properties of the disk 
(jZhu et al.l 120101 ) . The outbursts predicted in the IZhu et al.l (120101 ) study would have been 
most numerous during the very early stages (first few ~10 5 years of evolution), and decayed 
in intensity and frequency with time. The outbursts would have led to larger temperatures 
within the disk, allowing greater volumes of dust to be annealed than when more "normal" 
disk evolution was occurring. As outbursts were largely limited to early times, this would 
again suggest that the greatest amount of dust annealing would occur during this early 
epoch of disk evolution. The efficiency at which such materials are transported outward and 
preserved in the disk, however remains to be evaluated. The differences between the a=10~ 3 
and a=10 -4 runs in this study suggest that a modest decrease in the rate of disk evolution 
reduces the outward mobility of small dust particles. As such, the amount of crystalline 
silicates found in the cold, outer regions likely depends on the number of outbursts a given 
protoplanetary disk experienced, with larger numbers allowing greater amounts of dust to be 
processed and for greater levels of outward transport to occur. This, however, is somewhat 
speculative and should be the focus of future work. 

Another issue to explore in future studies is the vertical motions of the solid materials 
during transport and how they would affect the inferred mineralogy of disks, since obser- 
vationally, we are only able to see the very surfaces of a protoplanetary disk. An issue to 
consider is whether the dust at the disk surface layer would be representative of the bulk 
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mineralogy of the disk interior, where planet formation is expected to occur. Vertical diffu- 
sion w o uld car ry processed grains up to the surface of the disk where they could be observed. 
Cieslal (l2010bl ) showed that at a given location in the disk, for the case of a vertically uniform 
a, the dust in the upper layers of the became representative of the entire column of dust on 
a timescale of t ~ H 2 /D z where H was the local disk scale height and D z the local vertical 
diffusion coefficient. For standard simplifications of D z =ac s H , this gives t where Q is 

the local orbital frequency. The story becomes more complicated, however, if the vertical dif- 
fusivity varies with height above the disk midplane, as would be exp ected in disks where the 
magn eto- rotational instability (MRI) is either active throughout (e.g. lFromang fc Papaloizou 



1996 



2006) or if it is limited to just the surface layer over an other wise "Dead Zo ne" (e.g. iGammie 



Turner et al. 


2010; 


Zhu et al. 


2010) 



grains in upper layers would become "trapped" there by the turbulent eddies for long periods 
of time, spending extended periods of time in these high altitudes before returning back to 
the deeper regions of the disk. As such, there would be less direct communication between 
the disk midplane and the surface layers, offering the possibility of the chemical and min- 
eralogical species seen at high altitudes differing from that what was present in the larger 
fraction of the disk. Further, here we used the simple ID a-viscosity disk model in order 
to illustrate how the particle-tracking techniques could be applied in an evolving disk. This 
model assumes no significant variation in the radial transport of materials with height above 
the disk midplane. Previous studies had shown that in a p urely viscous disk, radial transport 



varied significantly with height above the disk midplane ( )Urpinlll984j : iTakeuchi Sz Linl 12002 



Keller fc Gailll2004t ICieslal 120071 . 120091 ). In these models, o utwar d mo vement is mo re efficient 
at the disk midplane than the surface layers, and ICieslal (120071 ) and I Cieslal ( 120091 ) explicitly 
showed how this could lead to dif ferences in the prope rties of materials at the disk surface 
versus that in the interior. While iFromang et all ( 120111 ) has found such ordered flows do not 
necessarily develop in fully MRI driven disks, the possibility of variations in radial transport 
with height do remain. This issue requires a detailed understanding of the large-scale flows 
that develop in the disk, and their variations with time, and should be the focus of future 
work. 

The results of this work also have important implications for the use of the meteoritic 
record to infer the types of processing that occurred in our solar nebula. Given the wide 
variety of paths that solids may take in the solar nebula, and thus the range of conditions 
that a swarm of "sibling" particles (those that exist at the same location in the disk at a 
given time) may have seen, it is necessary to collect data on a large number of particles in 
order to say anything conclusive about the origin of the group as a whole with any confidence. 
We must identify what properties are shared among a majority of primitive objects in the 
meteorites and which are limited to a small fraction. Those rare or extreme properties 
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may be attributed to those particles that take uncommon paths through the nebula. Thus 
we must take caution in inferring properties of the solar nebula or the origin of primitive 
materials that are based on observations of a single or few particles. 

The author is grateful for detailed comments and suggestions from Joe Nuth and an 
anonymous referee that led to substantial improvements in this paper. This work was sup- 
ported by NASA Grant NNX08AY47G awarded to F.J.C. 
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Fig. 1. — Comparison of the particle-tracking model described here (dashed lines) with the 
dye-tracking model typically used in transport studies (solid lines). Plotted are the surface 
densities of trace materials in a steady-state disk as described in the disk, at times 10 4 (black 
lines), 10 5 (light grey), and 10 6 (dark grey) years. The bumps and wiggles in the surface 
density distribution at 10 6 years in the particle-tracking model arises due to the fact that 
many particles are lost from the disk with time, leaving a small number behind to define the 
distribution. 
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Fig. 2. — As in Figure 1, except for the case of the evolving disk as described in the disk. 
Here, just like in Figure 1, the good agreement between the two models demonstrates the 
particle-tracking model accurately describes the collective dynamical behavior of the particles 
in a protoplanetary disk. 
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Fig. 3. — The radial locations of 3 particles from the simulation from the evolving disk 
described above. These particles originate at the same location in the disk and finish the 
simulation (at t=10 6 years) at roughly the same location in the disk ~2.56 AU. As can 
be seen, while their point of origin and final location are similar, each particle takes a 
different path through the disk, which would lead to exposure to and residence in different 
protoplanetary disk environments. 
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Fig. 4. — Left: The paths of the particles shown in Figure 3 during the first 10 5 years. 
Right:The corresponding temperatures seen by these particles in that time. 
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Fig. 5. — Plotted are the extreme radii (outermost =triangles and innermost=crosses) for 
each of the surviving particles in the evolving disk simulation described above plotted as a 
function of the final position of the particles at the end of the simulation. Left: All particles 
remaining in the disk. Right: Those particles whose final position ranges from 1 to 5 AU. 
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Fig. 6. — Maximum temperature seen by each surviving particle plotted versus the final 
position of those particles at the end of the simulation. Left: All particles remaining in the 
disk. Right: Those particles whose final position ranges from 1 to 5 AU. 




Fig. 7. — Same as Figure 5, except the turbulence parameter is assumed to have a value of 
a=l(T 4 . 
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Fig. 8. — Same as Figure 5, except the turbulence parameter is assumed to have a value of 
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Fig. 9. — Annealing timescales for the three annealing rates considere d here: Mg smokes from 
Fabian et al.l (120001 ) (FabOO: solid line), Mg glass as described by Bockel ee-Morvan et al. 
J2OO2J) (BM02: dashed line), and Mg smokes as studied by lHallenbeck et al.l (boooj ) 
(HalOO:dash-dot line). 
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Fig. 10. — Left: The fraction of each grain that became crystalline (annealed) plotted as 
a function of the peak temperature seen by the grains, using the pa t hs for those particles 
in the evolving disk described within the text and the iFabian et al.l ( 120001 ) annealing law. 
Right: The fraction of each grain that became crystalline plotted versus its final location in 
the disk at the end of the simulation. 
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Fig. 11. — Same as Figure 10, except using the iBockelee-Morvan et al.l (120021 ) annealing law 
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Fig. 1 2. — Same as Figure 11, except using the Silicate Evolution Index of lHallenbeck et al 

feoooh . 
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Fig. 13. — The crystallinity distribution of particles of radius 10 /im, 100 /im, and 1 mm, run 
with t he same disk model and a s sump tions described above, using the intermediate annealing 
law of iBockelee-Morvan et al.l (120021 ) . Larger grains, regardless of their crystallinity, are 
confined to radial locations closer to the central star due to the effects of gas drag (note the 
different scales on the x-axis). Despite the decrease in the survival frequency of the larger 
particles, the fraction of grains that are crystalline among the survivors is the same for all 
sizes considered here. 
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Fig. 14. — The time that each grain became crystalline plotted versus the final location o 
the grain for the three annealing laws considered here. Note that grains annealed late in 
disk history are limited to the inner disk, whereas those that make it to the outer disk come 
from a relatively small time window in the history of the disk. 



